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Abstract 

For the analysis of the r-mode oscillation of hot young neutron stars, it is 
necessary to consider the effect of differential rotation, because viscosity is not 
strong enough for differentially rotating young neutron stars to be lead to uni- 
formly rotating configurations on a very short time scale after their birth. In 
this paper, we have developed a numerical scheme to solve r-mode oscillations 
of differentially rotating polytropic inviscid stars. This is the extended version 
of the method which was applied to compute r-mode oscillations of uniformly 
rotating Newtonian polytropic stars. By using this new method, we have suc- 
ceeded in obtaining eigenvalues and eigenfunctions of r-mode oscillations of 
differentially rotating polytropic stars. 

Our numerical results show that as the degree of differential rotation is 
increased, it becomes more difficult to solve r-mode oscillations for slightly 
deformed configurations from sphere compared to solving r-mode oscillations 
of considerably deformed stars. One reason for it seems that for slightly 
deformed stars corotation points appear near the surface region if the degree 
of differential rotation is strong enough. This is similar to the situation that 
the perturbational approach of r-mode oscillations for slowly rotating stars in 
general relativity results in a singular eigenvalue problem. 



1 Introduction 

Recent intensive as well as extensive investigations of the r-mode instability have 
revealed its important role in neutron star physics (!§],]§]; for a review see e.g. [|[ ||; 
for r-mode oscillations see |& @, R). At the early stage of those investigations, 
the r-mode oscillations were found to become significantly unstable for some range 
of the core temperature ||, 0, [H^. However, models used in those papers were 
too simplified to deduce a decisive conclusion for real neutron stars because one 
could only obtain r-mode oscillations of uniformly and slowly rotating Newtonian 
polytropic stars and apply the results to realistic neutron stars. 

Actual neutron stars are more complicated configurations. There are many fac- 
tors to be included to get a definite answer about the effect of the r-mode instability 
on them. For example, neutron stars are general relativistic objects, they are gener- 
ally expected to be rapidly rotating at their birth, the physical states of the interior 
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cannot be quantitatively approximated by polytropic relations, the surface of the 
neutron stars is not considered to be fluid, they have magnetic fields and so on. 

The rotation law is one of such factors. Although in almost all papers about 
the r-mode instabilities neutron stars are assumed to rotate uniformly, newly born 
neutron stars whose ages are less than 1 year should be considered to be rotating 
differentially. Some neutron stars are newly formed from core collapses of massive 
stars on a dynamical time scale without enough time to redistribute the angular 
momentum due to viscosity and rotate differentially. (At some stages of recycled 
neutron stars, they may begin to rotate differentially due to instability of f-mode 
oscillations via gravitational wave emission because angular momentum increase due 
to accretion onto old neutron stars leads them to bifurcation points to configurations 
within which there exists internal motion. However, this instability results in non- 
axisymmetric configurations and they are out of the scope of the present paper. ) 
Thus we need to study the effect of differential rotation of newly born neutron stars 
on the r-mode oscillations. 

Moreover, some authors discuss the onset of differential rotation of a star due to 
self-induction by r-mode oscillations. Spruit |11 argued that the back-reaction of 
gravitational radiation by the r-mode oscillation on the stellar fluid may induce the 
differential rotation of the star. Rezzolla et al. |l2|, |1| 14 studied a drift induced 
by the r-mode oscillation of the stellar fluid which is initially uniformly rotating, 
and concluded that the r-mode instability may lead to the differential rotation of 
the star. By using numerical simulation of stellar hydrodynamics, Lindblom et 
al. |l5[ ] studied non-linear evolution of an initially uniformly rotating star by the r- 
mode instability and found that the stellar fluid develops strong differential rotation. 
Thus in the general context of the r-mode instability, it seems to be important to 
investigate the r-mode oscillations of differentially rotating stars. 

In our present investigation, by extending the numerical method which we devel- 
oped to compute r-mode oscillations of uniformly and rapidly rotating polytropes 
|l6| , |Tij |, we have succeeded in obtaining sequences of r-mode oscillations for dif- 
ferentially rotating polytropic stars (in a different context, the r-mode oscillations 
of configurations with slightly differential rotation are treated in |l^, [l^]). From 
our numerical results for differentially rotating stars, we have found that, if the 
degree of differential rotation is small, r-mode oscillations of the ordinary type can 
be obtained. Here r-modes of the ordinary type mean that the eigenfrequency has 
discrete spectrum. However, when the degree of differential rotation becomes suffi- 
ciently large, no r-mode oscillations of the ordinary type seem to be allowed, because 
in the surface region of the star corotation points appear and the basic equation 
becomes singular there. 

Since the purpose of this paper is to show the characteristic feature of the r- 
mode oscillation of differentially rotating stars, we only present a limited results of 
numerical computations, i.e. for the m = 2 mode, where m is the azimuthal mode 
number of the perturbation. 



2 Formulation of the problem 

We study the effect of differential rotation on the r-mode oscillations by solving 
the linearized equations of fluid motions of stars. The scheme employed in this 
paper is the extended version of that developed to analyze the r-mode oscillation 
of uniformly rotating Newtonian polytropes |T^, Therefore, we will explain the 
scheme only briefly. 
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2.1 Assumptions on unperturbed stars 



We construct axisymmetric equilibrium configurations of differentially rotating poly- 
tropic inviscid stars in Newtonian gravity by employing the Straight-Forward-Newton- 



Raphson method 20 . Here the polytropic relation is expressed by the following 
formula: 

p = Kp 1+ ^ 1 (1) 

where p, p, K and N are the pressure, the density, the polytropic constant and the 
polytropic index, respectively. In this paper polytropes with N = 0.5, 1.0 and 1.5 
are investigated. 

As for the rotation law of differentially rotating unperturbed stars, two kinds of 
angular velocity, ft, distribution are employed: 

n= ^ ( 2 ) 

and 

n= , M . (3) 

where R is the distance from the rotation axis, R eq is the equatorial radius of the 
star and f2 c is the central angular velocity. The quantity A is a parameter which 
represents the degree of differential rotation. The rotation becomes more differential 
as A becomes smaller. On the other hand, when we consider the limit of A — > oo, 
the stellar rotation tends to uniform rotation in both cases. For the region where 
R/R eq ^> A, the former law, Eq. (Q), tends to that of constant specific angular 
momentum distribution, while the latter, Eq. (|^), tends to that of constant linear 
velocity distribution. Thus we will call the former as a j-constant rotation law and 
the latter as a v-constant rotation law, hereafter. 



2.2 Linearized basic equations for perturbed stars 

As mentioned before, the scheme to handle the r-mode oscillation of differentially 
rotating polytropes is the extended version of that to solve the perturbed fluid equa- 
tions for uniformly rotating polytropes by the Newton-Raphson iteration scheme 

MM- 

Our basic equations consist of linearized parts of the following equations: 1) 
the continuity equation, 2) the ^-component of the vorticity conservation equation, 
or the compatibility equation between the r- and (^-components of the equation of 
motion, 3) the ^-component of the vorticity conservation equation, or the compati- 
bility equation between the r- and ^-components of the equation of motion, and 4) 
the (^-component of the equation of motion. The linearized gravitational potential 
in the integral form is used in the (^-component of the equation of motion. Here, 
the spherical coordinates (r, 6, (p) are used. 

Concerning the perturbed quantities, we assume the harmonic expansion of the 
physical quantities as 

Sf(r, 0, (p,t) = ^2 exp(-icrt + im<p)f m (r, 6), (4) 

m 

where / denotes a certain physical quantity, 5 is the Euler perturbation of the corre- 
sponding quantity, a is the frequency of the oscillation and m is the mode number. 
Because of the stationary and axisymmetric nature of unperturbed stars the equa- 
tions are separable with respect to the variables t and ip. Thus our basic equations 
are characterized by a and m. The perturbed fluid equations are expressed as 
follows: 
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where r s (9) is the surface radius of the equilibrium configuration and Sr s is the 
change of the surface. Here, functions f n (r,r') are defined as 



fn(r,r') 



(r>>r), 



(10) 



and pq denotes the mass density of the equilibrium states, whereas Su r , Svg, 
Sw v , Sp and 5<j) are the Euler perturbations of the r-, 9- and ^-components of the 
velocity in the corresponding orthonormal frame, the density and the gravitational 
potential, respectively. 

In these equations, we have made use of the following adiabatic perturbation 
between the pressure and the density: 



6p = ( l\5p 
P \ Nj p> 

where Sp is the Euler perturbation of the pressure. 



(11) 
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2.3 Boundary conditions 

Since we treat infinitesimal oscillations around equilibrium configurations, all the 
perturbed quantities must behave regularly throughout the whole space. In partic- 
ular, the Eulerian change of the density and the Eulerian change of the gravitational 
potential must be regular functions of the position. 

As for the perturbed flow velocity, the boundary condition on the surface can 
be expressed as 

§r ( 6ur - h^ SV8 ) + {mn a)6p = °- (12) 

This condition comes from a requirement that the fluid element on the equilibrium 
surface is displaced to the perturbed surface. 

3 Numerical results 

We can obtain eigenvalues and eigenfunctions by solving numerically the basic equa- 
tions mentioned above by the Newton-Raphson iteration scheme, together with the 
boundary condition. In actual computations, we use 44 mesh points in r-direction 
and 11 mesh points in ^-direction. 

In Figure 1, the eigenvalues of equilibrium sequences with different values of 
the parameter A are plotted against the ratio of the rotational energy, T, to the 
absolute value of the gravitational energy, W, for the m = 2 mode of TV = 1 
polytropes. In this figure, for the sake of simplicity, the eigenvalues are normalized 
by using the central angular velocity. This choice of the normalization constant 
results in decreasing of the eigenfrcquencics for the smaller value of A or for the 
larger degree of differential rotation. It is clearly seen that the eigenvalues decrease 
as the ratio T/\ W\ increases for a given value of A. Terminal points of sequences at 
the larger end of T/| W| correspond to final models beyond which our present code 
could not give converged solutions. 

Figure 2 shows ^-component of the perturbed velocity on the equatorial plane 
for the m = 2 oscillation of N = 1 polytropes with different degrees of nonuniform 
rotation. The horizontal axis is the normalized stellar radius. The equilibrium stars 
are rotating rather rapidly, i.e. r p = 0.7. Here r p is the axis ratio which is defined 
by 

r p = r ax /r cq , (13) 

where r ax and r cq are the stellar radius along the rotational axis and that on the 
equatorial plane, respectively. 

Three curves in Figure 2 correspond to equilibrium models of rigid rotation, 
differential rotation with the parameter A = 1.3 and extremely differential rotation 
with A = 0.6, respectively. The values of T/|W| are almost the same for these 
models, i.e. T/\W\ = 0.073,0.082 and 0.080 for the uniformly rotating model, the 
model with A — 1.3 and the model with A = 0.6, respectively. Thus we can regard 
these three models as those with the same rotational amount but different degrees 
of differential rotation. From this figure, we can see that if we take the differential 
rotation into consideration, oscillatory motions of the stellar fluid are confined only 
in the narrow layer near the surface of the star for a model with larger degree 
of differential rotation. In other words, for differentially rotating models, if the 
degree of differential rotation is large, only the fluid near the surface region can 
oscillate appreciably but the inner bulk of the star almost remains at its original 
position. This behavior can be seen not only in the equatorial plane but also in all 
9— directions. 

In Figure ^, eigenfrequencies of r-mode oscillations for the m = 2 mode of 
differentially rotating polytropes with several polytropic indices are plotted against 
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the value of T/|W|, i.e. N = 0.5,1.0 and 1.5. Here we have used the rotation 
law (|^) with A = 1.3. Terminal points of theses sequences are final models beyond 
which our code could not give solutions to the linearized basic equations because 
of very slow convergence. From this figure, it is clear that, although the values of 
a /£l c are not the same, the tendency of decrease of the normalized eigenfrequency 
as increase of T/|VF| is the same. 

For oscillations of differentially rotating stars, we need to pay attention to emer- 
gence of corotation points at some places in the stellar interior because appearance 
of corotation points corresponds to a singular behavior of the basic equations (see 
e.g. |2l[). Here corotation points are defined as points at positions where the 
following condition is satisfied: 

& = <T-mn = 0. (14) 

The condition of appearance of corotation points can be expressed by a simple 
inequality for differentially rotating stars. For the j-constant rotation law, Eq. (^), 
the quantity a is rewritten as 



" c (^ (R/Z^ + Ai) ■ (15) 



If this quantity a becomes zero within the star, we cannot solve the eigenvalue prob- 
lem because the matrix which governs the problem becomes singular as mentioned 
above. From Eq. (|l5|), the quantity a takes its minimum value on the rotational 
axis and its maximum value on the equatorial surface. Hence, 

^-m)<a-mn<n c (^--^). (16) 

For m = 2 r-modes of N = 1 polytropes, as seen from Fig. 1, the lower limit is 
always negative in the star. Thus the condition for which the quantity a becomes 
zero can be obtained by requiring that the upper limit of a must be positive inside 
the star. As a result, the condition that corotation points exist in the star can be 
written as 

n c * TTA2' (17) 

where m = 2 is substituted because we only discuss m — 2 oscillation modes in this 
paper. In order to find out whether the corotation points appear in the star or not, 
we need to check this condition for oscillations of differentially rotating stars with 
the j-constant rotation law. 

In Figure 4, values of eigenfrcqucncics which satisfy the equality of the above 
condition are plotted against the value of log A. We will call this curve as a critical 
curve. If the eigenfrequency of a certain star is located above the critical curve, a 
corotation point appears in the star. On the other hand, if the eigenfrequency is 
below the critical curve, no corotation point appears. In this figure, eigenfrequen- 
cies of two equilibrium sequences, sequences with r p = 0.95 and with r p = 0.70, 
are also plotted. As seen from this figure, the eigenvalues locate well below the 
critical curve when the degree of differential rotation is small. When we decrease 
the value of A, i.e. increase the degree of differential rotation, however, the eigen- 
frequencies approach to the critical curve. We can infer from the result here that, 
for extremely slowly rotating stars for which the standard slow-rotation approxima- 
tion is applied, the permitted range of the parameter A for the existence of discrete 
modes is narrower than that for rapidly rotating stars. 

We show the ratio of the eigenfrequency to the critical value for appearance of 
corotation points in Figure 5. From this figure, it is clear that eigenfrequencies are 



6 



approaching toward the critical curve monotonically as the value of A is decreased. 
For nearly spherical configurations, i.e. the sequence with r p = 0.95, the eigenfre- 
quency of the model with A — 0.615 reaches up to 97% of the critical value. On the 
other hand, for rapidly rotating models (r p = 0.70), eigenfrequency reaches only 
93% of its critical value, when A = 0.585. 

For slowly rotating stars (r p ~ 1), sooner or later, the eigenvalues seem to 
cross the critical curve, though it is very difficult to get converged solutions of 
oscillations by using our numerical code in such cases. This is not related to the 
ability of the code because r-mode oscillations can be easily obtained for nearly 
spherical configurations with larger values of A. 

For the v-constant rotation law, the condition that corotation points of m = 2 
mode appears in the star is written as 

->^=- (18) 

In Fig. 6, values of eigenfrequencies which satisfy the equality of the above con- 
dition are plotted against the value of log A as in Fig. 4 of the j-constant rotation 
law. The curves denoted by rapid rotation and slow rotation correspond to the 
eigenfrequencies of the equilibrium sequences with r p = 0.6, and 0.95, respectively. 

As seen from this figure, the eigenfrequency reaches only up to about 71 % of 
the critical frequency even for the model with the same axis ration of r p = 0.95, 
A = 0.25 and T/\W\ = 0.010. It should be noted that this value of T/\W\ = 0.010 is 
smaller than the value of T/|W| = 0.031 for the model with r p = 0.95 and A = 0.65 
for the j-constant rotation law. It implies that corotation points are more likely to 
appear for the rotation law with a steeper angular velocity distribution from the 
rotational axis to the equatorial surface. 



4 Discussion and summary 

4.1 Influence of differential rotation on the spin evolution of 
young hot neutron stars 

It is important to know how far differential rotation will change the spin evolution 
scenario of young hot neutron stars by the r-mode instability. In order to get 
precise answer to this problem, we need to perform time evolutionary computations 
of rotating stars. However, it is a difficult task to do even for uniformly rotating 
stars. For the estimation of spin evolution of uniformly rotating neutron stars, 
therefore, only the growth rates of the unstable modes have been investigated by 
evaluating the time scales of r-mode instability due to gravitational emission and 
damping times of the modes due to viscosity (see e.g. || [|). 

Although, for uniformly rotating stars, the scheme to evaluate time scales has 
been already established by using the energy of the mode as measured in the rotating 
frame p2| , p3| , which is defined as a quadratic form of Eulerian perturbations, 
there is no corresponding definition of the energy functional which can be used to 
determine the stability of the oscillation modes for differentially rotating stars. Thus 
in this paper we will discuss only a tendency of evolution of differentially rotating 
stars via instability of the r-mode. 

There are two main reasons why the r-mode instability for the uniformly rotat- 
ing stars becomes important for spin evolution of young hot neutron stars. First, 
the growth rate of the instability due to gravitational wave emission is so large that 
the spin change of neutron stars would occur on a very short time scale. Second, 
due to the characteristic feature of r-mode oscillations, the stabilizing effect due to 
(bulk) viscosity is much weaker than that for the f-mode oscillations. Although it 
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is not easy to estimate quantitatively the effect of these two mechanisms for differ- 
entially rotating stars, we note that the time scale of gravitational emission from 
differentially rotating stars may be larger compared with that for corresponding 
uniformly rotating stars. 

The energy loss rate due to gravitational wave emission can be written as 



dE _ 1 
~dt ' ~32tt ^ 

l>2,l>m 



{l+1) D, 



(19) 



where Di m and J; m are the mass multipole moments and the current multipole 
moments, respectively |m| . Here < -' +1 - ) -D/ m and J/ m are the (£ + l)-th time 
derivatives of the corresponding quantities. These two multipole moments are de- 
fined as: 

D lm = J r l pYrd 3 x , (20) 

2 1 
~cl + l 



Jim = -t^T / r l (pv)-(rxX7Yr)d 3 x , (21) 



(22) 

where c is the speed of light and Y™ are the spherical harmonics. Since the energy 
radiated away from the system is a large fraction of the rotational energy of stars 
(measured in the inertial frame), E Iot , the time scale of the gravitational wave 
emission, tqw, can be roughly estimated as 

This estimation can be applied to all rotating stars irrespective of their rotation 
laws. 

As seen from the eigenfunctions of the perturbed velocity fields, oscillations with 
large amplitudes for differentially rotating stars are confined to the surface region 
where the density is small. On the other hand, the cigcnfrcquencies of differentially 
rotating stars are not different significantly from those of uniformly rotating stars. 
Thus the time derivatives of the current multipole moments of differentially rotating 
stars are considerably smaller than those of uniformly rotating stars. It implies that 
the time scale tqw for differentially rotating stars would be longer than that esti- 
mated for uniformly rotating stars. Therefore, the instability of r-mode oscillations 
of differentially rotating stars would not be so significant as that of corresponding 
uniformly rotating stars having the same amount of rotational energy. Of course, 
this must be checked by computing the growth rates of oscillations quantitatively 
by using such a scheme as was employed to analyze the f-mode instability of rapidly 
rotating stars Jj5|. 



4.2 R-mode oscillation of differentially rotating stars and the 
strange nature of the r-mode oscillation of the perturba- 
tional approach in general relativity 

As discussed in Introduction, r-mode oscillations of neutron stars should be analyzed 
in the framework of general relativity. Although r-mode oscillations in general 
relativity have been investigated by several authors, there appear some features 
which do not exist in uniformly rotating Newtonian stars ]26| , 
of them is related to the strange nature of the eigenvalue problem if barotropic 
or isentropic configurations are treated. Some authors conclude the eigenvalue 
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problem becomes singular pq, p7fl a nd other authors claim that the system results 
in an over-determined state ]29|. . 

The strange nature of the eigenvalue problem for barotropic stars in general rel- 
ativity arises from appearance of the following term in the coefficient of the highest 
order derivative term of the master equation if the slow rotation approximation is 
employed: 

-on^-mfn- \ {e+ \> » ), (24) 

where u>(r, 6) is the dragging of the inertial frame of rotating general relativistic 
configurations and £ is the index of spherical harmonics. 
If we introduce the following effective angular velocity 

2(n-a;(r,0)) 

then we can write the quantity ctgr as follows: 

ctgr = <7 - mtt c s . (26) 

Here, even if the stellar angular velocity, O, is constant value, the effective angular 
velocity, Q e s depends on the location. Thus the peculiar nature of the eigenvalue 
problem in general relativistic slow rotation approximation can be essentially the 
same phenomena of appearance of corotation points for differentially rotating stars 
which has been discussed in the previous section. 

In Figure 7, distribution of the quantity Qeg/Q for £ = m = 2 on the equatorial 
plane is shown for uniformly rotating general relativistic polytropes with ./V = 1.0. 
Relativistic polytropes are defined by the following relation: 

p = Ks 1+1 ' N , (27) 

where e is the energy density. The strength of gravity is measured by the ratio 
of the pressure to the energy density at the center Q. Models shown in Fig. 7 
are configurations with N = 1 and p c / e c = 0.5 which corresponds to the typical 
parameter for neutron stars. Here p c and e c are the central values of the pressure 
and the energy density, respectively. Two curves correspond to equilibrium models 
with r p = 0.97 and 0.59. As seen from this figure, although rotation speed of two 
configurations is quite different, the quantity il e g/Q shows little difference. Conse- 
quently, if the eigenfrequencies of the r-mode for general relativistic configurations 
behave similarly as those in Newtonian models, in other words, if the eigenfrequen- 
cies tend to decrease as the stars rotate rapidly, almost the same behavior as those 
of differentially rotating Newtonian stars may be expected. 

Therefore, there arises a possibility that the quantity <tgr can be negative 
throughout the star and does not change the sign for rapidly rotating configu- 
rations as differentially rotating Newtonian stars do. In other words, the ordinary 
type r-mode oscillations may be possible for rapidly rotating general relativistic 
stars. 

Moreover, as our present numerical results show, if the corotation point is about 
to appear in the star, it becomes very difficult to solve ordinary type r-mode os- 
cillations, i.e. r-mode oscillations with discrete eigenfrequencies and this seems 
consistent with the general relativistic analysis under the slow rotation approxima- 
tion. Furthermore, the eigenfunctions for such configurations behave like a delta 
function as seen from Fig. 2. 

We may need to explain the reason why some authors |^9|, |3(| could solve r-mode 
oscillations for isentropic stars in the post-Newtonian approximation. As our nu- 
merical results show, weakly or intermediately differentially rotating configurations 
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do not suffer from appearance of corotation points. Thus for weakly relativistic 
configurations, there can be the ordinary type r-mode oscillations, although they 
do not have exactly the same character as that of Newtonian configurations [^9[ |3(j . 

From our "numerical experimental data" in this paper, there is a possibility of 
obtaining ordinary type r-mode oscillations for rapidly rotating general relativistic 
configurations, if one will be able to develop a scheme to handle linearized equations 
of the general relativistic fluid equations as well as the linearized Einstein equation 
for rapidly rotating stars. 

4.3 Summary 

We have succeeded in developing a new scheme which can be used to solve r-mode 
oscillations for stars which rotate differentially. This is a necessary development in 
order to study the scenario that neutron stars have lost their angular momentum 
via gravitational waves just after their birth by the r-mode instability. 

By analyzing the r-mode oscillations of differentially rotating polytropes, we 
have shown that the eigenvalue problem is likely to become singular for configura- 
tions with large degree of differential rotation. 
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Figure 1: Eigcnfrcqucncies for m ~ 2 mode oscillations of N = 1 poly tropic equi- 
librium sequences with different values of A are plotted against the ratio of the 
rotational energy to the absolute value of the gravitational energy, T/|W|. The 
attached number to each curve is the value of A along the sequence. The eigenfre- 
quency is normalized by using the central value of the angular velocity tt c . Terminal 
points of the curves for larger value of T/IIFl correspond to final models beyond 
which our numerical code could not give solutions of r-mode oscillations. 
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Figure 2: The ^-component of the perturbed velocity on the stellar equatorial plane 
is plotted against the distance from the center for N = 1 polytropes and the m = 2 
mode oscillation. Three curves show distributions of the ^-component of the fluid 
velocity with different degree of differential rotation: the solid curve for a rigidly 
rotating configuration, dashed curve for a configuration with A = 1.3 and the dotted 
curve for a configuration with A — 0.6. The values of the axis ratio for the three 
models are the same, r p = 0.7, and the values of T/\W\ are T/\W\ = 0.073,0.082 
and 0.080 for the uniformly rotating model, the model with A = 1.3 and the model 
with A = 0.6, respectively. It should be noted that the large amplitude oscillations 
are confined to a narrow surface region for the highly differentially rotating model. 




Figure 3: Same as Figure 1 but for different values of N. Attached values to the 
curves are the polytropic indices, N = 0.5, 1.0 and 1.5. The rotation law is that of 
the j-constant law with A = 1.3. 
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Figure 4: The critical curve for appearance of corotation points and eigenfrcquencies 
of two equilibrium sequences are plotted against the value of log A Solid curve 
denotes the critical curve. Dashed and dotted curves show the eigenfrequencies for 
the equilibrium sequences with r p = 0.95 (slow rotation) and 0.70 (rapid rotation), 
respectively. Terminal points of eigenvalue curves correspond to models with A = 
0.615 for the r p = 0.95 sequence, and A = 0.585 for the r p = 0.70 sequence, 
respectively. The values of T/|W| along the sequences are not exactly but roughly 
the same. For r p = 0.95 sequence, T/\W\ = 0.023 - 0.031 and for r p = 0.70 
sequence, T/|VF| = 0.078 <~ 0.080. These terminal points of the curves for smaller 
values of A correspond to final models beyond which our numerical code could not 
give solutions of r-mode oscillations. 
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Figure 5: The ratio of the eigenfrequency to the critical value for the appearance 
of corotation points is plotted against the value of A The maximum value of this 
ratio reaches 97% of the critical value for the appearance of corotation points along 
the r p = 0.95 sequence. 
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Figure 6: Same as Fig. [| but for the v-constant rotation law. Dashed and dotted 
curves show the eigenfrequencies for the equilibrium sequences with r p = 0.95 and 
0.60, respectively. The values of T/|W| along the sequences are T/|W| = 0.010 ~ 
0.011 for the r p = 0.95 sequence and T/\W\ = 0.109 - 0.113 for the r p = 0.60 
sequence. 
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Figure 7: Distribution of the quantity fl e g/fl is shown on the equatorial plane of 
uniformly rotating general relativistic polytropes with N = 1.0. The strength of 
gravity parameter is large, i.e. p c /s c — 0.5. The solid and dashed curves show the 
models with r p = 0.97 (slow) and 0.59 (rapid), respectively. 
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